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| The internal structure of the jets produced in neutral current interactions for 

Q 2 > 125 GeV 2 has been studied using the subjet multiplicity with the ZEUS 
detector at HERA. Jets are identified in the laboratory frame by applying the 
^ | \ longitudinally invariant fey-cluster algorithm. Next-to-leading order QCD cal- 

^vq I culations have been obtained and compared to the data; a good agreement 

04 ■ has been observed in the region where hadron-to-parton corrections are small, 

I . ET.jet > 25 GeV. In this region, the strong coupling constant is extracted and 

^ [ the resulting value is a s (Mz) = 0.1185±0.0016 (stat.)^o'oo48 ( s y st -) -0 0071 (th.). 

■ 

o : , . 

t— 1 , 1 Introduction 

The study of the internal structure of jets gives insight into the transition between 
partons produced in a hard process and the experimentally observable spray of 
hadrons. This type of analysis is usually made using jet shapesu, where the energy 
flow inside the jet is considered, or the subjet multiplicitycl, where the study is 
done with jet-like structures (subjets) within a given jet. These structures are 
called subjets and are formally defined as the original jets, but resolved at a smaller 
I scale. At sufficiently high jet energy, where fragmentation effects are small, the 

internal structure of the jets is expected to depend mainly on the parton-radiation 
^ ■ pattern and should be calculable in pertubative QCD. 

Studies of the internal structure of jets have been made in pp collisions at 
Tevatroni and in e + e~ interactions at LEP@. It has been observed that phenomeno- 
logical models for the parton radiation based on QCD give a good description of 
the data. In addition, the results are in agreement with the expectation that the 
internal structure is determined by the type of the primary parton. It is observed 
that gluon-initiated jets are broader (i.e. contain a larger number of subjets) than 
quark-initiated jets, as predicted by QCD from the larger color charge of the gluon. 
At HERA, measurements of subjet multiplicities have Dreviously been presented in 
quasi-real photon proton collisions (photoproduction)H and in neutral current deep 
inelastic scatteringQ. In photoproduction, the mean subjet multiplicity (< >) 
is observed to become larger as the jet pseudorapidity (r]j et ) increases, in agree- 
ment with the predicted increase in the fraction of gluon-initiated jets. In deep 
inelastic scattering for jets defined in the Breit frame, it has been observed that 
jets pcontain a smaller number of subjets with increasing transverse energies and 
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towards the forward region of the Breit frame. 

New measurements of the mean subjet multiplicity for jets produced in neutral 
current deep inelastic scattering with the ZEUS detector at HERA are presented 
here. The jets are selected using the longitudinally invariant fcy-cluster algorithmi 
in the laboratory frame. Next-to-leading order (NLO) QCD calculations for < 
n s bj > are now available for jets defined in this frame and they are compared to 
the data after hadronization corrections are applied. From this comparison, the 
value of the strong coupling constant is extracted. 

2 Measurement of the mean subjet multiplicity 

The data sample used in this analysis was collected with the ZEUS detector at 
HERA and corresponds to an integrated luminosity of 38 pb _1 . During 1996- 
1997, HERA operated with positrons of energy E e = 27.5 GeV colliding with 
protons of energy E p = 820 GeV. Neutral current events were selected by requiring 
a positron candidate identified in the main calorimeter (UCAL). After removing 
the background contributions, events with Q 2 > 125 GeV 2 were kept and the 
jet algorithm was applied to the UCAL cells belonging to the hadronic system. 
The jet search is performed in the pseudorapidity (77) - azimuth (0) plane of 
the laboratory frame. The jet variables are defined according to the Snowmass 
convention!. The inclusive sample of jets with transverse energy Exjet > 15 GeV 
and — 1 < r]j et < 2 has been studied. 

For each jet in the sample, the jet algorithm was re-applied to all the parti- 
cles assigned to the given jet. The clustering was stopped when the distances yij 
between all pair of particles i, j are above some resolution scale y cut 

min\ET i, Et j} 2 / A 9 A , 9 \ 

to = — 1 J£" ( Ar] i + A(f) i) > v**- 

^T,jet 

All clusters found in this way are called subjets. 

The quantity ^ycutErjet determines the scale at which the internal structure 
of the jet is resolved. If this quantity is large, then the subjet structure is given 
by parton-radiation processes in which the relative transverse momentum is hard 
enough and perturbative QCD should be applicable to describe the subjet structure 
of the jets. 

The mean subjet multiplicity is defined as the average of the number of subjets 
in the considered sample of jets at a given y cut . The mean subjet multiplicity has 
been studied as a function of y cut and of the jet variables for a fixed y cut = 0.01. 

A simulation of the detailed properties of the hadronic final state is available in 
various Monte Carlo event generators. The generated samples of events have been 
passed through a detector simulation and analysed by the same program chain 
as the data. The response of the detector to jets of hadrons and the correction 



a The pseudorapidity is defined as 7/ = — ln(tan((9/2)) where 9 is the polar angle with respect to 
the proton beam direction. 
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Figure 1. Measured mean subjet multiplicity corrected to hadron level for jets produced in neutral 
current DIS interactions as a function (left) of the jet pseudorapidity and (right) of the jet 
transverse energy in the laboratory frame. 



factors for < n s bj > have been determined by applying the same jet algorithm to 
the hadron level of the simulated events. 

Figure |1] shows the measurements of the mean subjet multiplicity as a function 
of the jet pseudorapidity and jet transverse energy for y cut = 0.01. The data have 
been corrected to the hadron level and compared to the predictions of two Monte 
Carlo models. The color-dipole model as implemented in ARIADNEa gives a very 
good description pi the data while the matrix element plus parton shower (MEPS) 
model of LEPTO&3 is above the data. 

It is observed that the mean subjet multiplicity has no rjj et dependence. This 
is different than the results observed in photoproductioni, where an increase is 
observed as r)j et increases. These two different behaviors can be understood in 
terms of the admixtures of quark- and gluon-initiated jets in the two reactions. In 
the analysis of neutral current deep inelastic scattering presented here, the observed 
distribution is consistent with the expectation that the sample is dominated by 
quark-initiated jets from the Born process for the entire r]j et range. Comparisons 
with the predictions by the Monte Carlo models for quark- and gluon-initiated jets 
confirm this explanation. 

On the other hand, the mean subjet multiplicity is observed to decrease as the 
transverse energy of the jets increases. This means that the jets become narrower 
as Exjet increases, in agreement with the expectations. A similar behavior has 
been observed in previous results at HERA, Tevatron and LEP by using subjet 
multiplicities or jet shapes. 



3 NLO QCD calculations 

For jets defined in the laboratory frame, it is possible to obtain predictions for the 
mean subjet multiplicity from perturbative QCD to the second order in a s . The 
perturbative QCD calculations have been obtained with the program DISENTH. 
The number of flavors was set to 5 and the renormalization (/j,r) and factoriza- 
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Figure 2. Measured mean subjet multiplicity corrected to hadron level for jets produced in neutral 
current DIS interactions. The data are compared to NLO QCD calculations for different values 
of a s : (left plot) < n s bj > as a function of y C uu (right plot) < n s t,j > at y cut = 0.01 as a function 
of the transverse energy of the jets in the laboratory frame. 

tion (/xp) scales were both set to Q. The main source of theoretical uncertainty 
is due to terms beyond NLO and has been estimated by varying between Q/2 
and 2Q, keeping \i F fixed at Q. 

To compare with the data, the predictions given by DISENT (at parton level) 
have been corrected for hadronization effects by using the Monte Carlo models. In 
the ARIADNE and LEPTO models, the hadronization is performed by using the 
LUND string; mo dell as implemented in JETSElM 

The perturbative QCD predictions corrected for hadronization effects are com- 
pared to the data in figure 0. The predictions give a good description of the 
data, even where parton-to-hadron corrections are large. NLO QCD predictions 
have been obtained for different values of a s . The parametrizations of the proton 
parton distribution functions assuming different values of a s {Mz) as given in the 
CTEQ4A-seriesE-3 have been used; this procedure takes into account the correla- 
tions between the value of a s [Mz) used in the matrix elements and that assumed 
in the parametrization of the parton distribution functions. 

The comparison of the data with the theoretical predictions for different values 
of a s shows that the measurements are sensitive to a s (Mz). 

4 Determination of a s 

The value of a s has been extracted in the region where parton-to-hadron correc- 
tions are below 15%. The method consists of making a parametrization of the 
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dependence on a s of the theoretical predictions for the mean subjet multiplicity. 
A preliminary study of the systematic uncertainties has been performed; the most 
important contribution is the uncertainty on the modelling of the hadronic final 
state. 

The extracted value of a s (Mz) using < n s ^ > at y cut = 0.01 and for jets with 
E Ttjet > 25 GeV is 

a s (M z ) = 0.1185 ±0.0016 (stat.)^E (syst.)^;S (th.) 
in good agreement with the world average valued. 



References 



1. S.D. Ellis, Z. Kunszt and D.E. Soper Phys. Rev. Lett. 69, 3615 (1992). 

2. S. Catani et al, Nucl. Phys. B 377, 445 (1992) and Nucl. Phys. B 383, 419 
(1992); M.H. Seymour, Nucl. Phys. B 421, 545 (1994) and Phys. Lett. B 
378, 279 (1996). 

3. CDF Collab., F. Abe et al., Phys. Rev. Lett. 70, 713 (1993). 
DO Collab., S. Abachi et al., Phys. Lett. B 357, 500 (1995). 

DO Collab., V.M . Abazov et al., preprint FERMILAB-PUB-01-248-E (2001), 
|hep-ex/0 1080541 . 

4. ALEPH Collab., D. Buskulic et al., Phys. Lett. B 346, 389 (1995). 

5. ZEUS Collab. Contributed paper to ICHEP'99 (Tampere, July 1999) N-530. 

6. HI Collab., C. Adloff et al., Nucl. Phys. B 545, 3 (1999). 

7. S. Catani et al., Nucl. Phys. B 406, 187 (1993); S.D. Ellis and D.E. Soper, 
Phys. Rev. D 48, 3160 (1993). 

8. J. Huth et al. Proc. of the 1990 DPF Summer Study on High Energy Physics, 
Snowmass, Colorado, edited by E.L. Berger (World Scientific, Singapore, 1992) 
p. 134. 

9. L. Lonnblad, Corap. Phys. Comm. 71, 15 (1992); version used: 4.08. 

10. G. Ingelman, A. Edin and J. Rathsman, Comp. Phys. Comm. 101, 108 
(1997); version used: 6.5. 

11. S. Catani and M.H. Seymour, Nucl. Phys. B 485, 291 (1997). 

12. B. Andersson et al, Phys. Rep. 97, 31 (1983). 

13. T. Sjostrand, Comp. Phys. Comm. 39, 347 (1986); T. Sjostrand and M. 
Bengtsson, ibid. 43, 367 (1987). 

14. H.L. Lai et al., Phys. Rev. D 55, 1280 (1997). 

15. Particle Data Group, D.E. Groom et al., Eur. Phys. J. C 15, 1 (2000); S. 
Bethke, J. Phys. G 26, R27 (2000) 



5 



